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ABSTRACT 

We examine neutron star properties based on a model of dense matter com- 
posed of B=l skyrmions immersed in a mesonic mean field background. The 
model realizes spontaneous chiral symmetry breaking non-linearly and incor- 
porates scale-breaking of QCD through a dilaton VEV that also affects the 
mean fields. Quartic self-interactions among the vector mesons are introduced 
on grounds of naturalness in the corresponding effective field theory. Within a 
plausible range of the quartic couplings, the model generates neutron star masses 
and radii that are consistent with a preponderance of observational constraints, 
including recent ones that point to the existence of relatively massive neutron 
stars M ~ 1.7 'Mq and radii R ~12-14 km. If the existence of neutron stars 
with such dimensions is confirmed, matter at supra-nuclear density is stiffer than 
extrapolations of most microscopic models suggest. 

Subject headings: neutron stars, equation of state, skyrme model 



1. Introduction 



Ne utron star astronom y initiated by the serendipitous discovery of the first radio 
pulsar ( Hewish et al. 19681), ha s since found over 1700 similar "rotation-powered" neutron 
stars (lATNF Pulsar Catalogue! ) . Pulsar timing measurements in radio binaries yield a simple 
(unweighted) mean neutron star mass (M) ~1.4M Q . In contrast, measurement of general- 
relativistic parameters in neutron star-white dwarf binaries suggest (M) ~1.6M , owing to 
a few exceptionally large inferred neutro n star masses in t he latter case. Examples include 



the binary component PSR J0621+1002 flNice et al. 20071 ) with 1.7Z i 6 M Q at the la level 
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and a more stringent value of M >1.68M (p at the 2a level set by a neutron star binary in 
the Terzan 5 cluster (jRansom et al. 20051 ). Constraining neutron star radii is harder, due to 
uncertainties in atmospheric modelling and distance estimates, but bounds from the rmally 
emitting neutron stars eg. RXJ 1856.5-3754 (IWalter fe Lattimer 20021 ; iHo et al. 20071 ) imply 
that the canonical range of 10-12km is exceeded. These observations have oriented attention 
towards 'atypical' neutron stars with large mass and possibly large radius. 

An independent determination of mass (2.1O±O.28M ) and radius (13.8±1.8km) of the 
bursting ne utron star in the low-mass X-ray binary (LMXB) EXO 0748-676 has been claimed 
(jOzel 20061 ). b ased on an accurate determination of gravitationally red-shifted Fe and O ab- 
sorption lines (jCottam et al. 2002l ). However, the evolving nature of the source has compli- 
cated further observational tests of the same. Furthermore, Doppler tomography of emission 
lines in the mass-transfer stream between the neutron star a nd its less massive c ompanion 
suggests a more canonical value of 1.35M© for the former (IPearson et al. 20061 ) . NASA's 
upcoming Constellation-X mission will improve on such spectral measurements, shrinking 
systematic errors on mass and radius even further. With due caution on the observational 
front, the confirmation of a mass ~2.OM for a neutron star strongly constrains the equa- 
tion of state of dense matter, ruling out the possibility of extreme softening at high den- 
sities. It also implies an upper bo und on the energy density of observable cold and dense 
matter (ILattimer fe Prakash 20051 ). A large radius R ~13-14 km for a 1.4M© neutron star 
implies a stiff symmetry energy at densities of [1-2] no, where hq is the saturation density of 
nuclear matter. These connections between the nuclear physics of dense matter and neutron 



star observations have been the focus of recent reviews (jSteiner et al. 20051 ; ISedrakian 2006 
Page fe Reddv 20061 : Eattimer fe Prakash 2007h . 



Uncertainties in theoretical aspects of many-body interactions at n>[l-2]n lead to pre- 
dictions for the mass versus radius curve that vary widely depending on the equation of state 
(EoS), with the ma ximum mass M max ranging from 1.4-2. 7M & and radius at maximum mass 
-Rmax from 9-14km (ILattimer fe Prakash 20011 ). Constraints fro m astrophysical obse rvations 
and terrestrial laboratory data have whittled this range down f ILi fe Steiner 20061 ) so that 
most microscopic models of longstanding for nuclei and nuclear matter would struggle to 
explain the existence of relatively heavy neutron stars (2M ) which also have a large ra - 
dius ( -R~13km). In fact, on ly few stiff equations of state eg. , MPAl (IMuther et al. 19871 ). 
MS0 jMiiller fe Serot 1996h and PALI JPrakash et al. 1988ah are consistent with the EXO 
0748-676 constraint at the la level, and even they fail if the matter accreted onto the neutron 
star is helium-rich 0. However, while satisfying astrophysical constraints, it is also important 



1 The largest source of systematic error in extracting the mass and radius of EXO 0748-676 comes from 
the accreted mass fraction of Hydrogen 0.3 < X < 0.7 (Ozel 2006). 
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to keep in mind constraints from laboratory data on strongly interacting matter around nu- 
clear saturation density. This point was nicely brought out in recent papers by Li & Steiner 
(2006), and Klahn et al. (2006). 

In this conte xt, our goal in this paper is to explore further a r ecently proposed model of 



a skyrmion fluid (jOuyed fc Butler 19991 ; iJaikumar fc Ouyed 20061 ) , henceforth referred to as 



OBJ, that was shown to lead to a very stiff equation of state and consequently generate a large 
maximum mass as well as radius for a neutron star. However, it was pointed out (Lattimer, 
private communication) that the rapid rise of the compressibility and symmetry energy just 
above saturation de nsity in this model puts it at odds with experimental constraints from 



collective flow data flDanielewicz et al. 20021) and isospin diffusion studies in medium-energy 
heavy- ion collisions (ITsang et al. 20041 ). In this work, we determine the extent to which 
we can satisfy these constraints by extending the skyrmion fluid model to include higher- 
order interactions among the vector mesons that are theoretically motivated by arguments 
of naturalness in the corresponding effective field theory. The mass and radius predictions 
of the extended model, henceforth referred to as JOM, are also confronted with constraints 
set by the observation of X-ray burst oscillations, kiloHertz quasi-periodic oscillations in 
LMXBs and thermal emission from neutron stars. We include observational uncertainties 
wherever they may impact our conclusions. Our phenomenological model is able to satisfy a 
preponderance of these constraints. We emphasize at the outset that our model, in its current 
form, has not been investigated for its applicability to nuclei or more complex phases of 
matter at sub-saturation densities. Our EoS presently applies only to infinite nuclear matter 
and neutron-rich matter in the range [1-5] no- 

This paper is presented as follows. In § [21 we discuss some conventional EsoS for 
neutron stars; in §[3] we revisit the Ouyed-Butler-Jaikumar (OBJ) model for skyrmion stars 
and motivate higher-order interactions that serve to tune the stiffness of the equation of 
state such that laboratory constraints are met. The main features of the mass versus radius 
curves are explained in § HI We compare the results obtained in the skyrmion star model 
to predictions of other neutron star models in light of observational bounds in § [51 Our 
conclusions are in § [61 



2. Equations of State for Neutron Stars 

An equation of state for dense matter is a relation between pressure and energy (or 
baryon) density, usually derived from an underlying microscopic model or effective the- 
ory for strong interactions. To apply to neutron stars, it should be able to generate at 
least a 1.4M static neutron star with a radius in the 10-14km range. The connection 
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to the underlying microscopic theory can be formulated in several ways: examples include 
relativistic mean field th eory, non-relativistic po t ential models, relativistic D irac-Brueckner- 
Hartree-Fock theory etc ( lArnett fc Bowers 19771 ; lLattimer &: Prakash 20011 ) . As we are con- 
cerne d with recent finding s of rel atively heavy neut ron stars, we consider three s tiff eos: 
MSO jMiiller fc Serot 1996h . APR jAkmal et al 1998h and UU Jwiringa et al. 198sh . 



(i) Miiller & Serot (1996) used a relativistic theory of point-like nucleons interacting 
via mesonic degrees of freedom. These are the neutral scalar (a) and vector (u) fields, plus 
the isovector p meson. In this model, like or unlike meson-meson interactions are encoded 
by terms that are polynomials in the fields. By demanding a match to the properties of 
nuclear matter at saturation, they obtained a sequence of EsoS that depend on the coupling 
constants of the polynomial interactions. The stiffest eos (MSO) corresponds to vanishing 
couplings and yields a maximum mass of 2.7 M & . This model is consistent with a large 
neutron star mass and radius ~14km for static configurations. Presently, the la limits on 
the radius of the bursting neutron star source EXO 0748-676 are 13.8±1.8km. However, there 
is no fundamental symmetry principle that requires the higher-order couplings to vanish. 
Introducing natural values for these couplings drastically reduces the maximum mass to 
< 2.OM . 

(ii) Ak mal et al. (1998) obtai ned the APR EoS based on the Argonn e v-\ s nucleon-nucleon 
interaction (IWiringa et al. 19951). Urbana IX three - nucleon interaction (IPudliner et al. 19951 ) 
and a relativistic boost term (IForest et al. 19951 ) as microscopic input. This EoS gives a 
maximum neutron star mass of 2.2M and does not lie within the la limits on the mass- 
radius estimate of the neutron star in EXO 0748-676, even assuming the accreted matter is 
mostly hydrogen. 

(iii) EoS UU is obtained via similar variational methods applied to an older two-nucleon 
and three-nucleon interac tion (Urbana fu+U VII). This model yields a maximum mass of 
2.2 M & at radius 10 km (IWiringa et al. 19881 ). However, a radius larger than about 11km 
is not supported by this equation of state for a static neutron star with mass greater than 
1.4M . 

Our purpose in selecting and highlighting these EsoS is two-fold. Firstly, these models 
are representative of complementary philosophies behind constructing an equation of state 
for dense matter: as in (i) forgo the connection to laboratory data on vacuum two-nucleon 
interactions and focus instead on the empirical properties of large nuclei and infinite nuclear 
matter within a relativistic theory; or as in (ii) and (iii) insist on a satisfactory description 
of available data on the structure and interaction of few nucleon systems (free or bound) 
in a non-covariant approach. Secondly, these are among the stiffest equations of state that 
arise from hadronic degrees of freedom alone. It is possible that hybrid equations of state 
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that allow for quark matter at high density can be almost just as stiff (lAlford et al. 20071 ). 
but we will not consider quark matter EsoS in this work. 

These three EsoS are used in our neutron star mass-radius plots (FigsJUSJ) and compared 
with the EsoS for dense skyrmion matter which we now describe. 



3. The Skyrmion fluid 

3.1. Nuclear Matter Phenomenology 

Before the advent of Quantum Chromodynamics (QCD), T. H. R. Skyrme proposed a 
description of baryons as topological solitons in a meso nic field theor y that realizes spon- 



taneous chiral symmetry breaking in non- linear fashion ( jSkyrme 196ll ). This model is now 
qualitatively supported by studies of large N c QCD which suggest that mesonic degrees of 
freedoms are fundamental and baryons arise as solitons. When augmented by the inclusion 
of low-lying vector mesons (my <lGeV) and flavor symmetry breaking effects, the Skyrme 
model can provide a reasonable descriptio n of static baryon prope rties such as mass split- 



tings, charge radii and magnetic moments (ISchecter fc Weigel 20001 ) . In the 2-nucleon sector 



the problem of finding a sufficiently attractive isoscalar central and spin-orbit force within 
the Skyrme model at distances lfm< r <2fm has held up progress. One solution is to in- 
clude a dilaton field that mocks up scale-breaking in QCD and provides attraction in these 
channels. 

To make progress towards an equation of state for a skyrmion fluid, Kalbermann (1997) 
introduced a self-consistent model that incorporates medium effects through the response 
of the dilator]^ a and the isoscalar u-field to a smooth density distribution obtained by 
integrating over the collective co-ordinates of the Skyrmion. This is equivalent to an ensemble 
of non-interacting B=l Skyrmions, a valid picture upto a separation of 0.8fm (n <5no). A 



subsequent work (jJaikumar fc Ouyed 20061 ) extended this model to asymmetric matter by 
incorporating the p meson in the spirit of standard mean field approaches. The Lagrangian 
for the Skyrme model, augmented by the a, u fields, and including isospin-breaking effects 
from the p meson as well as explicit scale-breaking effects from the dilaton and quark masses 
is given in Jaikumar & Ouyed (2006). A fit to nuclear matter phenom enology is achieved 



through additional parameters in the dilaton potential, which is given by (jKalbermann 19971 ) 



2 Since chiral symmetry is broken non-linearlythe usual a field as the chiral partner of the pion does not 
appear in the theory. 
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V(a) = J B[l + e 4CT (4(7-l) + a 1 (e- CT -l) (1) 
+a 2 (e CT - 1) + a 3 (e 2CT - 1) + a 4 (e 3a - 1)] , 

where B « (240 MeV) 4 is related to the Bag constant (the non-perturbative glue that breaks 
scale-invariance in QCD). The six unknown parameters of the model are a\-a±, g w the u- 
N coupling and g p , the p-N coupling. To determine the a^'s, the following constraints are 
imposed: the scale anomaly condition ^| CT=0 = which implies that a\ = a 2 + 2a 3 + 3a4, 
the stationarity w.r.t ctq, viz., dE/dao = 0, a binding energy/nucleon of -16 MeV for infinite 
nuclear matter at saturation density (no = 0.16fm~ 3 ), and a choice of the compressibility 
K. E is the energy density of the fluid and o"o is the non-vanishing mean field value of 
the time component of the a field. At saturation, o"o is determined by the choice of the 
effective mass M, which then also fixes g w . The choice of symmetry energy and effective 
mass at saturation fixes g p . Once the a^s are determined, <To is generally obtained from 
its equation of motion for an arbitrary density. The a^s show very weak dependence on 
the choice of K in the range 200-300 MeV, while displaying more sensitivity to the choice 
of effective mass. Without further modifications, the model displays a sharp rise in the 
compressibility just above saturation, rising from K=240 MeV (our choice) to K ~ 2000 
MeV for a 10% increase in baryon density. Similarly, the symmetry energy rises too steeply in 
this range to be consistent with experimental constraints (see §3.41) . These inconsistencies are 
a consequence of the specific form of the dilaton potential, which is essential to preserve the 
trace anomaly relation (scale-breaking). Therefore, a modification of V(a) is not desirable. 
It is the exponential sensitivity of the curvature of the potential V{a) to the dilaton VEV 
that drives the compressibility to large values. One way to address this issue is to view the 
Skyrme model in the mean field approximation as an effective field theory of hadrons, so that 
the Lagrangian can be extended to include higher-order terms (meson self-interactions) that 
parameterize unknown physics at a more microscopic level. This rationale is also employed 
in Miiller & Serot (1996) although their model has an explicit a meson, point-like nucleons 
and no scalar-vector mixing, while our model has a dilaton as the only scalar, and includes 
scalar-vector mixing. If the meson fields are viewed as relativistic functionals, the higher- 
order interactions can be thought of as parts of an effective potential that determines their 
mean field values at a particular den sity through the stati onarity of the effective action 
associated to the Skyrme lagrangian (IFurnstahl et al. 19961 ) . They therefore modify the 
density dependence of the meson fields as well as the properties of the background skyrmion 
fluid that couples to these fields. We restrict ourselves to quartic self-interactions in the p 
and u fields with coupling constants whose value can be surmised by naturalness. Then, 
higher-order terms such as six or eight-meson self-interactions do not substantially change 
the results obtained in the quartic case. In the next section, we implement this procedure 
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to obtain an acceptable behaviour of the compressibility and symmetry energy in dense 
Skyrmion matter. 

3.2. Mean-field equations 

The additional quartic interactions take the form 



c u = £°+c nt ; C nt = +|^ 4 , (2) 

£° = —F^^ + ^mW-^mj (3) 

for the u) field and 

C p = C° p + £?- t C? = +% g y, (4) 

C° p = —F^ + ^my-^njp (5) 

for the p field. Here, the F's are the standard field-strength tensors for the abelian (lu) and 
non-abelian cases (p). Working at T = 0, we have 



(h.3 , h.3 \ ( h.3 _ h.3 

y^Fp + K F n J y K F p K F r 



where fc^ is the Fermi momentum of species i = neutron, proton. Then, the Skyrme La- 
grangian is compactly expressed as 

C = C 2 + £ 4 + C u + C p -V{a), (7) 

where £ 2 , £ 4 involve gradients of the Skyrmion profile. For the densities of interest (n < n < 
5no) where the approximation of non-overlapping skyrmions is valid, this profile drops off fast 
enough that the mean-field averaging simply counts the number of individual Skyrmions in 
a given volume, equivalent to a non-interacting Fermi gas model. At n > 5no, we expect cor- 
rections to our mean field model from Skyrmion overlap. We also do not expect such a mean 
field treatment to apply at densities much lower than saturation density, since Skyrmions do 
not form a uniform fluid there. Thus, our model is restricted to (n < n < 5n ). For the ui 
and p fields, the equations of motion read as follows: 

a^ul + b^ujQ + c u = 0; (8) 
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6 



,bu> 



0; 



a pPo + b pPo + c P 

4 

_ X9p , 2crn™ 2 



6 



e 2a °m;, c p 



-g w n . 



g p nS 



(9) 
(10) 

(11) 



where wn,Po denote mean field values and S—(l-2x), with x=(rii-n)/2n being the proton 
fraction of neutron-rich matter {x—1/2 for symmetric matter). The magnitudes § of the 
quartic couplings £, \ are estimated from the naturalness argument for an effective field the- 
ory, viz., that the co-efficients of the various terms in the Lagrangian, through a given order 
of truncation, should be of the same size when expressed in an appropriate dimensionless 
form. Thus, we find 



12 



e a °m u 

Mn 



192 



Mn 



(12) 



We choose the effective mass, given by M = e CT0 Mn, to be 600 MeV at saturation density, 
so that e <7 °=2/3 for a bare nucleon mass Mo =900 MeV (neglecting the ~40 MeV contribution 
from explicit symmetry breaking). Since the model is fit to saturation properties, g u itself 
depends on £, whose value must be chosen so as to satisfy the naturalness condition above. 
Furthermore, real solutions to the equation of motion for the dilaton Eqn. (1161) cease to exist 
beyond a small range of couplings. This restricts us to 0.1 < £ < 0.3 and 1.0 < x < 2.0. 
These values differ from those in Miiller & Serot (1996) due to additional factors of e 2cr ° from 
the dilaton (metric) that appear in the fitting expressions for g u and g p and the qualitatively 
different form of the dilaton potential (it contains all powers in a). The energy densities 
corresponding to the vector mean fields are 



4 2 
so the total energy density is then 



bppl 



c P Po 



E = E kin -+ 

Ekm = 



n,p 



E^ + E p + V{a) , 
'k F E F (E 2 F + k F ) 



(13) 



(14) 



M 4 



In 



k F + E F 
M 



where the kinetic energy E^ikp) comes from a Lorentz boost of the static Skyrmion to 
momentum kp (IKalbermann 19971 ). In Eqn.( fT4l) . E Fi = \ kp. + M 2 since the Dirac effective 



3 We choose the sign of the couplings to be positive since this guarantees zero mean fields at vanishing 
source density (baryon/isospin). 
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mass is the same for both neutrons and protons. The binding energy is E/n — M and the 
pressure is given by 



P = n 



0.5 ( E f I + 9^0 

V n,p J 



E. 



(15) 



The effective mass M is determined at any density from the equation of motion for o"o 



n,p 



kpEp 

4vr 2 



. ^ dV ,9,9 

+ 4E kiD + - b u uj 2 - h p pl = 

da 



(16) 



3.3. Compressibility 

For symmetric matter, 5 = and p vanishes. The compressibility is defined as 



an 



n 2 *Wn) + 2n d{E/n) 



dn 2 



dn 



In the present case, this is equivalent to 



(17) 



K = 9n 



d 2 E _ f d 2 E \ da 
V dnda J dn 



dn 2 

da _ ( d 2 E\ (d 2 E 
~dn ~ \dnda ) ' [ikr 2 

d 2 E 



+ 



9l 



dn 2 3nE F u u t o%^ 
d 2 E _ _2g u b u u 

Ep 
3nM 2 



dnda 
d 2 E 



da 2 



Ei 



-4 



dV d 2 V , 2 
2 + 6t> w u; 



O?0"O d<7j 







(6 W + g^gg) 



(18) 



(19) 
(20) 



For the choice K=240 MeV, binding energy=-16 MeV and M=600 MeV at saturation 
density, the best fit Oj values for various £ are listed in Table 1, correcting an unfortunate 
error in their values quoted for {=0 in Jaikumar & Ouyed (2006). 
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For £=0, at saturation, freedom in choosing the fit parameters a, allows for a deli- 
cate cancellation between various contributions to the compressibility such that K is small 
(~200-300 MeV). This is achieved largely as a result of fine-tuning d 2 V/da 2 away from its 
natural scale ~ B. As we move to slightly higher density, this fine tuning cannot be re- 
covered since V(ao) is exponentially sensitive to changes in the a VEV. Consequently, the 
contribution from the u;-meson appearing in the d 2 E/dn 2 term dominates, pushing the com- 
pressibility to unnaturally large values, as in the OBJ model. For £ 7^ 0, a w 7^ and u 
being large, it is clear from Eqn. fll9p that d 2 E/dn 2 and hence K is substantially reduced, 
as in the JOM model. This effect is reflected in the EoS of symmetric matter, shown in 
the upper panel of FigJTJ By definition, the slope of the pressure-density curve is pro- 
portional to the compressibility. The hatched region represents constraints on the EoS of 
symmetric matter from the analysis of coll ective flow data in nucleus-nucleus collisions at 
(l-2)GeV/nucleon (IDanielewicz et al. 20021 ) . While the £ = curve (OBJ) is clearly too 
stiff, choosing natural values of the coupling (0.1 < £ < 0.3) provides a considerable im- 
provement, and the JOM model is able to satisfy the flow constraint for £ ~ 0.25. The APR 
EoS, which has a phase transition at n ~2n , is also shown for comparison. 
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Flow constraint 
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Fig. 1. — Upper panel: EoS of symmetric matter for various values of the quartic coupling 
£ in the sk yrmion fluid model. The APR EoS is also shown. The hatched region is the flow 
constraint (IDanielewicz et al. 2002 ) . Lower panel: Symmetry energy 02 for natural values 
of the couplings £ and x m the same model. 
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3.4. Symmetry energy 



In neutron-rich matter, 5^0 and the energy/particle e = E/n can be expanded about 
the symmetric point 5 = (x = 1/2) : 



e(n, 5) = e{n, 0) + a 2 (n, a (n))5 2 + a 4 (n, a (n))5 4 + ... 
where ... are higher order terms, and 



(21) 



a 2 (n,a Q (n)) 



1 (d 



2! \dS 2 
1 /<9 4 e 



9 2 p P 3 f 



5=0 



127r 2 m2 e 2<J o 1 6£ F 



+ 



k 6 



648£ f 



4 3 3 
/cp kpEp Ep^ 



n" 



X 



3847 



(22) 



Studies of neutron matter show that the value of x for beta-equili brated matter obtaine d 
by retaining only a 2 is a good approximation over a wide range of n (IPrakash et al. 1988bl ). 
Recently, Steiner (2006) has examined the role of the quartic coefficient a 4 at high density, 
finding the effects on the EoS to be small, although the threshold for the onset of rapid 
neutrino cooling via the direct urea process can change considerably. As we are focusing 
on the equation of state in this work, we drop the a 4 term, keeping in mind that were we 
to retain the term in Eqn. (l21j) . the correction to the energy difference of symmetric and 
neutron-rich matter is at the 5-10% level for densities of interest. Proceeding with a 2 alone, 
as defined in Eqn. (j2"Tl) . beta equilibrium and charge neutrality conditions yield 



(3n 2 nx 



1/3 



4a 2 5 



(23) 



which fixes a solution x = xo(cr ) for any baryon density. Inverting this relation, and solving 
Eqn. fr23|) along with the mean field equations for a , u and p , we determine x (n) explicitly. 
For the density range n/n = [1..5], 0.045 < x$ < 0.11 for (,x / and the direct urea 
threshold is just about reached at the upper limit. This conclusion could be affected by the 
inclusion of the a 4 term in the expansion of Eqn. (|2ip . 

The lower panel of Fig{T]shows the density-dependence of the symmetry energy a 2 (o"o(n), n). 
Since a 2 in Eqn. (123!) depends on a , which takes different values for the same density in sym- 
metric and asymmetric matter, the curves do not begin at 32 MeV (our choice for symmetric 
matter). We have chosen a 2 from Eqn. fl23p to represent the symmetry energy since it is this 
quantity that determines the proton fraction of beta-equilibrated matter. This is different 
from the usual identification of the symmetry energy, which is made at 5 = 0. In addi- 
tion, it must be noted that a and hence a 2 now depend on x as well as £• F° r £=X=0, a 
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stiff symmetry energy results, that behaves approximately as 02 (?^) =02(^0) (^/no) L26 in the 
range [l-1.5]n . Li & Steiner (2006) have argued that such a parametric ally stiff symmetry 



energy is inconsistent with isospin-diffusion data i n heavy-ion collisions (jRami et al. 2000l ) 



and the measured neutron-skin thickness of Lead (IStarodubsky & Hintz 19941 ). For natural 
values of the couplings, the symmetry energy softens considerably, behaving approximately 
as a2(n)=a 2 (n ) (n/n ) ' 71 , which is consistent with the aforementioned experimental stud- 
ies. While increasing £ does not change the density dependence, increasing x softens the 
symmetry energy at high density. This is a consequence of scalar- vector mixing in our model. 

The pressure of neutron- rich matter in /^-equilibrium (excluding lepton pressure), scaled 
to that of a relativistic fluid with the same energy density, is shown in Fig. [2j Close to 
saturation, the pressure is well-described by P(n,5) = P Bym (n) + n 2 5 2 da 2 /dn where P sym is 
the pressure of symmetric matter. The curves in the lower panel of Fig [T] imply that 8 is 
smaller while da2/dn is larger for vanishing couplings, as compared to the case with non-zero 
couplings. This drives the initial rapid rise of the pressure relative to the energy density for 
£=X=0. For non-zero couplings, the softer symmetry energy and larger value of 8 lead to 
a larger pressure at saturation but a more gradual increase relative to the energy density. 
Increasing £ at fixed x softens the equation of state and has a progressively decreasing effect 
at high densities, where the kinetic contribution to the pressure begins to dominate; hence 
the gradual approach to the relativistic limit. Increasing x & t fixed £ has a more dramatic 
effect on the pressure and energy density, causing the EoS to exceed the relativistic limit at 
high densities. This is due to the dynamics dictated by the dilaton potential at high density 
and scalar-vector mixing in our model. The nucleon's effective mass increases rapidly with 
increasing density beyond a certain value of x an d the skyrmion shrinks, akin to a strong 
repulsive force. 
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Fig. 2. — Pressure of beta-equilibrated matter (relative to a relativistic gas of the same 
energy density) as a function of density for the OBJ and JOM EsoS. 
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Mass-Radius relation 



We employ the RNS code (IStergioulas fc Freidman 19951 ; IStergioulas fc Freidman 19981 ) 
to generate mass versus radius curves for a sequence of static and rapidly rotating neu- 
tron stars with the OBJ and JOM EsoS. In both cases, the composition of the star is 
as follows: (i) l<n/n Q <5: (n,p,e~) matter in beta-equilibrium with nucleonic pressure 
and e nergy density given by Eqns. (jT4l) and (|T5l) ; (ii) n/n <1: the BBP equation of 
state (IBaym et al. 1971al) for densities b elow nuclear saturation density, matched to the 



BPS equation of state (IBaym et al 1971bt ) for the low-density nuclear crust of the star. For 



non-zero couplings (JOM EoS), the maximum mass is not reached until densities larger than 
5n , so we had to extend our model to higher densities. Therefore, the maximum masses for 
JOM EoS indicated in Figj3] are to be viewed as extrapolations into a regime where the ap- 
proximation of non-overlapping B=l spherically symmetric skyrmions is not likely to hold. 
A more accurate method would have to first determine the topology and size of deformed or 
overlapping skyrmions, which we do not attempt in this work. 
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Fig. 3. — Mass-radius curves for static (left panel) and rotating (right panel) stars with the 
OBJ and JOM EsoS. 
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The variation in the mass-radius curves with the values of £ and x reflects the following 
two facts: the stiffer the EoS at supra-nuclear densities, the larger the maximum mass and 
(ii) the larger the pressure in the range [l-2]no, the larger the radius f or a (1-2) star. These 



correlations have been established and emphasized in previous work (jLattimer &: Prakash 20011 ) 
The extreme stiffness of the symmetry energy for £=x=0 results in large pressures in the 
range [1.5-2.0]no and leads to a large radius R ~ 14km for a 1.4M Q star in the OBJ model. 
The softening induced by non-zero couplings lowers the pressure in this range, leading to 
smaller radii in the JOM model. The maximum mass is also lower with respect to the case 
when the couplings are set to zero. For the rapidly rotating models, we chose a rotation fre- 
quency v ~600Hz, corresponding to a period of 1.6ms. In general, the additional centrifugal 
forces in a rotating star help to counteract the pull of gravity, resulting in larger radii for a 
given mass. 

As pointed out in the context of a different mean field model (Miiller & Serot 1996), the 
inclusion of the additional quartic terms softens the equation of state for beta-equilibrated 
neutron-rich matter considerably, making it difficult to obtain a neutron star mass larger than 
2M & . Within our model, while a large radius and mass is possible if these quartic terms 
are omitted, the requirement of respecting laboratory constraints near saturation density 
and arguments of naturalness imply that the inclusion of these couplings is essential and 
its consequences (a lowering of maximum mass and radius) quite general. We now turn to 
compare our results (along with those from other EsoS) with some current observational 
bounds on the mass and radius of neutron stars. 



-18- 



2.5 



1.5 



0.5 




10 



11 



12 13 14 

Radius (km) 



15 



16 



Fig. 4. — Observati onal constraints on neutr on stars: The region allowed by observations 
of RXJ 1856.5-375 4 ([Walter fc Lattimer 20021) is shown in red; regions excluded by glitche s 
in the Vela pulsar Link et al. 19991 ). bv thespin-rate of J1748-2446ad jrlessels et al . 200fih . 
and by causality are demarcated in pink. Mass-radius curves for static stars corresponding 
to some stiff equations of state (APR, UU, MSO) are plotted along with the ones in this 
work (OBJ, JOM). The z=0.35 constraint for EXO 0748-676 rules out the APR and UU 
equations of state at the la level but allows them at the 2a leve l. Yellow dashed lin es are 
90% confidence limits from observations of qLMXB X7 in 47Tuc (IHeinke et al. 20061 ). 
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Neutron star observations 



In Figs J4l5l we show several constraints on a neutron star's mass and radius that follow 
from general theoretical principles and observations of neutron star phenomena. Fig. H] is 
relevant for static or slowly rotating sta rs. A measurement of the blackbody radiation ra- 
dius -R 00 =16. 5km (ITrumper et al. 20041 ) a nd distance estimate rf= 117pc for the thermally 
emitting neutron star RXJ 1856.5-3754 (IWalter fc Lattimer 20021 ) yields the red allowed 
region, which is a lower bound on the radius since the blackbody radiator has the smallest 
emitting area for a given luminosity. Observations of glitches in the spin-down of the Vela 
pulsar place a one-parameter c onstraint on the fractional moment of inertia in the crust 
AI/I m 0.014 (ILink et al. 19991 ) . This parameter is the pressure at the core-crust interface, 
which takes typical values in the range 0.25 < P/(MeV/fm 3 ) < 0.65. The upper and lower 
limits translate to excluded regions in the upper left region of FigJD marked Vela-glitch 1 
and 2 respectively. Note that the region marked Vela-glitch 2, if taken at face value, rules 
out even a 1.6M star for the APR and UU EoS. The region marked Vela-glitch 1 still 
allows the APR and UU EoS up to 1.8M Q . Causality of the EoS excludes a (partially over- 
lapping) smaller pink region in the upper lef t corner. A scaling relatio n between mass and 
radius for the minimum allowed spin period (ILattimer fc Prakash 2004T). comb ined with the 
highest observed spin-frequency (716Hz) of J1748-2446ad (IHessels et al. 20061 ). excludes the 
region marked "rotation". 90% confidence limits on the radiation radius for the thermal 



source X7 in the globular cluster 47Tuc (IHeinke et al. 20061 ). arising from atmospheric mod- 
elling that includ es surface gravity effects consistently, are shown by the dashed yellow lines. 
Ozel's estimate (lOzel 20061) of mass and radius based on a red-shift measurement z=0.35 
(jCottam et al. 20021 ) and an assumed atmosphere of accreted hydrogen for EXO 0748-676 
gives the solid black line, with la and 2a limits displayed. 

We have plotted the mass-radius curve from the JOM and OBJ EsoS as well as the 
APR, UU and MS0 EsoS. Among these, only MS0 and OBJ lie well inside the red allowed 
region for a star~ 1.4M Q . A smaller distance, still within the above ±12 uncertainty lim- 
its, would also allow the JOM EoS, provided the mass of RXJ 1856.5-3754 is ~ 1.8M . 
The relatively soft equations of state such as APR and UU do not satisfy this constraint 
coming from RXJ 1856.5-3754 unless the mass of this object > 2.0M Q and the glitch con- 
straint is ignored. Although MS0 and OBJ appear promising in this light, as argued before, 
they are theoretically incomplete without higher-order terms, from the standpoint of an 
effective field theory. The addition of extra quartic terms to the Lagrangian in the OBJ 
model, which is essential for the self-consistency of the approximations and the truncation 
scheme, softens the EoS considerably. This yields the JOM EoS, which comes closest to 
satisfying all constraints. It is noteworthy that the APR, UU, OBJ and JOM mass-radius 
curves are also consistent with recently determined bounds (not shown in the figure) on 



-20- 



the radiatio n radius of neutron st ars in the globular cluster M13 (IGendre et al. 2003al ) and 
cu-Centauri fjGendre et al. 2003bl ) while MSO is not. Also, in addition to the JOM EoS, the 
APR and UU EoS are consistent at the 2a level, but not at the la level, with observational 
constraints set by X-ray bursts in EXO 0748-646 . 



Therefore, accepting all these observations as accurate, the OBJ and JOM EsoS are 
unique among the stiff EsoS considered here, in that they are most likely to satisfy all 
constraints from the afore-mentioned observations. Unlike the OBJ EoS however, the JOM 
EoS can also satisfy constraints from laboratory experiments, as demonstrated in §21 The 
JOM EoS also implies that the mass of the neutron star RXJ 1856.5-3754 has to be ~ 1.8M 
and that of the qLMXB X7 has to be > 1.6M Q . 
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In Figini we display mass-radius curves for rapidly rotating stars with the same EsoS 
as in the previous figure. Constr aints from kHz quasi-periodic oscillat ions (QPOs) in the 
bursting neutron star 4U 1728-34 (jMiller et al. 19981 ; iBarret et al. 20061) and practical upper 
limits on spin-frequencies of neutron stars with hadronic EsoS (ILattimer &: Prakash 2004 ) 
yield the wedge-shaped allowed region in blue. We have assumed a fairly typical QPO 
spin-frequency of 300 Hz in obtaining this constraint. We also employ compactness con- 
straints obtained from the analysis of Nath et al. (2002) which assume a two-spot model 
to explain the X-ray burst oscillations of LMXB 4U 1636-53. Their best-fit value for com- 
pactnes s is M/R=0.126 ass uming a spin- frequency of 290 Hz. With an assumed mass of 
1.6M Q (ICasares et al. 20061 ) for this object, we obtain the smallest of the allowed triangular 
region bounded by dashes in the central right region of the figure. With 90% (M/R— 0.163) 
and 99% (M/i?=0.183) confidence-level constraints, we obtain the larger triangular regions 
(also demarcated in Fig IS] by dashed lines). We plot the mass-radius curves for rotating con- 
figurations with z/=300Hz in order to make a fair comparison to the constraints. We observe 
that the rotating APR and UU configurations cannot satisfy the combination of constraints 
if the best-fit compactness value is used but the MS0, OBJ and JOM EsoS can. However, if 
we choose a lower mass M=1.44M Q or higher confidence limits for the compactness of 4U 
1636-53, the APR and UU EoS c an also satisfy thi s constraint. Accretion argues against 
having a mass smaller than 1.6M (IGiles et al. 20021 ) . Thus, an accurate mass measurement 
of this object, although complicated due to its accreting nature, would be very desirable. 

We have not attempted here to exhaustively survey all the EsoS that have been inves- 
tigated in the literature and compared to observations, nor utilize all possible observational 
con straints. There are a variety of EsoS based on re lativistic extensions o f potential mod- 
els ( jSteiner et al. 20051 ). simple parameterized EsoS (Prakas h et al. 1988aT) and relativistic 
field theories with density dependent couplings (IKlahn et al. 20061 ). only a few of which 
(eg. PALI, MPAl mentioned in the introduction) can achieve consistency with a majority 
of o bservational an d experimental data. There are constraints from neutron star seismol- 
ogy (IDuncan 19981 ) which are satisfied only with the softest EoS possible in our model, one 
that cannot generate a 2.OM maximum mass for a neutron star. While improved lim- 
its from future observations would naturally prove more definitive, the above discussion is 
meant to demonstrate that, at present, only an equation of state that is quite stiff at high 
density and moderately so near saturation density can satisfy a preponderance of astrophys- 
ical and laboratory constraints. The JOM EoS, based on a relativisitic mean field theory 
with scale-breaking and symmetry breaking inspired from QCD, provides one such promising 
example. 
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Fig. 5. — Constraints for rapidly rotatin g stars: The blue wedge is the regi on allowed 
by kHz QPO oscillations in 4U 1728-34 jMiller et al. 19981 : barret et al. 2006h assuming 
a typical spin frequency of 300 Hz, while the dashed wedges are the regions allowed by 
compactness values M/R=0.12Q (best- fit), 0.163 (90%) a nd 0.183 (99%) inferred from X-ray 
burst oscillations of LMXB 4U 1636-53 jNath et al. 20021 ). We assume its mass to be 1.6M Q . 
All mass-radius curves are for stars rotating with spin freqency u=300 Hz (period P=3.3ms). 
See text for details. 
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6. Conclusions 

We have examined a mean field model for dense nuclear and neutron-rich matter, with 
direct application to neutron star interiors. The model is derived from the Skyrme La- 
grangian, with the inclusion of the dilaton VEV performing a dual role: (i) to mock up 
scale-breaking in QCD and (ii) to dial the interaction between the density-dependent vector- 
meson mean fields and the skyrmion fluid which makes up the dense medium, thus achieving 
self-consistency. We included quartic self-interactions for the vector mesons to correct an 
unnaturally rapid rise of the compressibility and symmetry energy at densities just above 
saturation. The resulting equation of state is in better agreement with laboratory constraints 
on the compressibility and symmetry energy. At densities much above saturation, the equa- 
tion of state is stiff; consequently it yields a large maximum mass for a neutron star. This 
is particularly interesting given accumulating data that points to the existence of neutron 
stars with M ~ 2M , though corresponding systematic errors tend to be larger than is the 
case for stars with M ~ 1.4M Q . 

The pressure and symmetry energy of dense matter, as well as the mass-radius curves for 
neutron stars depend sensitively on the (natural) values of the quartic couplings. Increasing 
the u;-meson self-coupling reduces the compressibility as well as the symmetry energy near 
saturation density considerably, while increasing the p-meson self-coupling stiffens the EoS of 
neutron-rich matter at high density, thereby increasing the maximum mass. The maximum 
neutron star mass in our model lies between 1.8 to 2.0M Q for static configurations with 
corresponding radius at maximum mass between 10. 5-11. 5km. For rapidly rotating stars, 
these values increase to 2.O-2.3M and 11.5-12.5 km respectively. The minimum spin period 
for the maximum mass star lies between 0.66 ms and 0.71 ms. These macroscopic effects are 
traceable to changes in the dilaton VEV, which effectively controls the Skyrmion size through 
its interactions with the density-dependent vector-meson mean fields, thereby determining 
the stiffness of the EoS, and highlighting the role of scale-breaking in our model. 

Confronting our relatively stiff EoS with a set of observational constraints on neu- 
tron star mass and radius, we find encouraging agreement. Including commonly used EsoS 
based on extrapolations of microscopic models of the nucleon-nucleon interaction, we infer 
that matter is likely to be stiffer at supranuclear densities than such models would sug- 
gest. Therefore, an accurate and simultaneous determination of mass and radius for the 
more 'extreme' neutron stars, viz., those presently suggestive of high mass and large radius, 
will be particularly valuable. In addition, constraints from collective flow data in heavy-ion 
collisions and the isospin dependence of the strong interaction near saturation density are 
proving to be valuable benchmarks for dense matter studies. 
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Table 1. Fit parameters for the JOM model for various £ 7^ 0. The OBJ model has £=0. 
In both cases, fsT=240MeV, B.E =-16MeV, E sym =32MeV, g p =4M7 



£ 9uj ai a 2 a 3 a 4 



0.0 7.54 -1.699 -25.780 29.237 -11.464 

0.1 6.67 -1.960 -27.051 29.798 -11.502 

0.2 6.05 -2.260 -28.700 30.674 -11.636 

0.3 5.57 -2.602 -30.747 31.878 -11.870 



